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ABSTRACT 

The variability of 30 PG quasars has been observed in the red band during three 
years. A rest-frame structure function analysis shows an increase of variability with 
the time interval up to 0.2 mag r.m.s. after two years. A comparison with the 
lUE data available for PG quasars shows that the variability increases with the 
rest-frame frequency at each time interval. Evidence for this effect is supported by 
the analysis of the average variability of other published quasar samples. The effect 
can completely account for the increase of variability with redshift found in other 
studies. 

Subject headings: Galaxies: Quasars: General 
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1. INTRODUCTION 

The study of the quasar variabihty is of primary import tool for con- 

straining the physics of the emission (Rees 1984, Terlevich et aL 1992, Hawkins 
1993), and as a selection criterion (Trevese et al. 1989, Veron & Hawkins 1995). 
In both cases, it is essential to understand on statistical grounds the dependence of 
variability on the intrinsic luminosity of the source, on the redshift and the wave- 
length. For practical reasons, most variability studies with short time sampling are 
concentrated on the optically violent variables (OVVs), whose luminosity changes by 
a large amount over short periods of time. However, they represent a small fraction 
of the entire population. The majority of quasars selected in flux limited samples 
exhibit, in the optical band, variations of a few tenths of magnitude, with typical 
time scales of months to several years (see e.g. Usher, Warnock, Green 1983; Pica et 
al. 1988). In these samples selection effects can be more easily taken into account. 
Thus, we have observed a subsample of the Palomar-Green (PG) quasars (Schmidt 
& Green 1983), which are selected by a uniform and objective ultraviolet excess cri- 
terion. Our subsample of 30 objects has been observed with the Schmidt telescope 
60/90 cm of the Astronomical Observatory of Rome, equipped with a CCD detector 
(Natali & Pedichini 1990), which provides a field of 17' X 25', where several bright 
objects can be used for comparison. Photometry has been obtained in the R band, 
to optimize the efficiency, with sampling intervals ranging from about one month to 
3 years. A statistical comparison with archive lUE data has also been performed 
to explore the dependence of the variability on the wavelength. Such a dependence 
has been investigated by Cutri et al. (1985) who directly compared UBVRIJHK 
observations of a sample of 7 quasars, finding some evidence of a hardening of the 
spectrum during the bright phase. Such a behaviour has also been found by Kinney 
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et al. 1991 from lUE observation of a subsample of PG quasars. Giallongo, Trevese 
and Vagnetti (1991, paper I) have suggested that this behaviour may be responsible 
for the correlation of the amplitude of variability with redshift, found in their anal- 
ysis. In this paper we provide a new statistical evidence for the increase of quasar 
variability with the rest-frame frequency. 

2. OBSERVATIONS AND DATA ANALYSIS 

The 60/90 cm Schmidt telescope of the Astronomical Observatory of Rome is 
located at Campo Imperatore, on the Apennine mountains, at 2200 m above sea 
level. The focal length is 180 cm and the pixel size of the 384x576 Thompson 7882, 
UV coated, CCD chip is 23 x 23 /im. This provides a scale of 2.64 arcsec/pixel. 

Since the limiting magnitude of the PG sample is ~ 16.5, in each CCD field of 
about 17' X 25' centered on a QSO there are 20-40 stars within one magnitude from 
the QSO, which also appears star-like. This allows a good frame-to-frame calibration 
and a precise evaluation of the photometric accuracy. Short individual exposures 
of about 10 minutes were adopted for the brightest quasars to avoid saturation. 
The resulting total exposure times range from 10 to 20 minutes. Each object was 
observed three times on average. The sample is presented in Table 1, marked as 
CIMP (Campo Imperatore). Twilight fiat fields were taken for each observing run. 

A catalog of reference objects in each quasar field was built up using a simple 
automatic detection algorithm. The object positions were refined computing the 
baricenter of the light distribution. At each epoch, relative R magnitudes were 
computed for three different apertures with radii ranging from 3 to 7 pixels. 

Eor each target QSO a multiple catalog was then derived, which contains the 
photometry at all the epochs in the three apertures. Magnitude offsets respect to 
a reference epoch, which was usually the oldest, were computed excluding the QSO 
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itself and those objects with magnitude difference > 3(7, which are less than fO% 
in the worst case. The average noise, associated with each pair of frames, has been 
defined as the standard deviation (Tij of the magnitude differences of the reference 
objects in the frames i and j. The optimal aperture for each pair of frames was 
selected as the one giving the lowest noise level. The fraction of light lost does not 
affect the relative photometry, since it is the same for all the stellar images, including 
the quasars. The distribution of noise estimates in the optimal apertures is shown 
in Fig. 1. The modal value, which is representative of the noise standard deviation 
(see Kendall & Stewart 1963, eq. 11.41), is ctr ~ 0.035. The resulting intrinsic 
variations 5m of the QSO magnitude were computed for each target and for each 
time interval Atrest = At /[I + z) (see below) and are shown in Fig. 2 and reported 
in Table 2, where columns 1 and 2 indicate the a and 5 coordinates, column 3 the 
redshift, column 4 the B magnitude from Schmidt and Green (1993), column 5 the 
observing time as day, month, year, column 6 the magnitude variation respect to 
the reference epoch, and column 7 the corresponding standard deviation (Tij. 



Rather than studying in detail the light curves of a few individual objects, it is 
possible to derive some general constraints on the emission mechanisms from the 
ensemble properties of statistical samples. The ensemble statistics is described by 
the "structure function" which is defined in different ways in the literature. Here 
we define the structure function as 



3. VARIABILITY ANALYSIS 



a) Optical Data 
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where the angular brackets indicate the ensemble average, At is the time interval 
between two observations, the exponent p may take different values, and Cp is a 
suitable normalization factor introduced to measure the structure function in units 
of standard deviation. For p = 2 and Cp = 1, eq. (1) corresponds to the first-order 
structure function of Simonetti, Cordes & Heeschen (1985). In this case, when the 
magnitude variations are due to intrinsic variability plus noise, the two contributions 
are added in quadrature, and the variability structure function becomes 5*2 = S2 — cr^, 
where cr^ is the noise variance. For p = 1, the average is more robust, because it is 
less sensitive to atypical data. Under the assumption that the statistical distribu- 
tions of noise and variability have the same shape, the factor Ci has the same value 



for both, with Ci = y7r/2 in the case of a gaussian distribution. We then adopt in 
the following analysis the variability structure function defined as: 



A similar analysis can be performed using the median instead of the ensemble aver- 
age, as done by Hook et al. 1994. 

If the structure function is computed as a function of the observed time intervals, 
it depends on the redshift distribution of the quasar sample. Thus it is physically 
more significant to compute the rest frame structure function after reduction of 
the time intervals to the rest-frames of individual quasars Atr = At/(1 -|- z). The 
structure function is shown in Fig. 3. The noise (T„ is computed as the rms value of 
the standard deviations (Tij associated with the frame pairs z,j contributing to the 
relevant bins. Fach error bar represents the statistical uncertainty on the average 
of the magnitude differences contributing to the same bin. 

At the short time scales ^100 days the variability is only slightly greater than 
noise, thus it is only marginally significant at the level of accuracy allowed by our 
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photometry. However, a systematic increase corresponding to average variations of 
about 0.1 mag after 1 year is clearly present. 

b) lUE Data 

To compare our observations of variability in the R band with other wave- 
lengths, we have examined the lUE archive through the European Space Infor- 
mation System (ESIS). We have extracted the available short-wavelength and long- 
wavelength spectra of the PG objects, excluding: (i) objects classified as Seyfert 
galaxies {Mb > —23.5), (ii) spectra with low S/N 1, per pixel), (iii) single epoch 
observations. The resulting sample consists of the 21 objects marked in Table 1. 
Eor about half of these objects, variability in the red is also available. Since the ob- 
servations of 3C 273 are much more numerous than those of other objects, we have 
used a randomly extracted subset of observations with comparable time sampling. 
Erom both the short and long wavelength spectra we have computed broad band 
magnitudes in the ranges 1300 — 1900 A and 2500 — 3000 A, with an arbitrary zero 
point. Erom these data we have computed the relevant structure functions (eq. 2), 
reported in Eig. 3. It appears that at each time interval the amplitude of variability 
increases with the rest-frame frequency, indicating, on average, a hardening of the 
spectra in the brighter phases, at least under the assumption that the power-law 
shape of the spectra is preserved. It should be noted that the latter hypothesis would 
imply that the structure functions at different wavelengths differ by a multiplicative 
factor. The presence of different time scales, increasing with wavelength, cannot be 
excluded by the present data. Notice that a better knowledge of the shape of the 
structure function in different wavelengths would allow to constrain possible devi- 
ations from power-law spectra, even in the absence of simultaneous multifrequency 
observations. 
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Our results are consistent with the previous suggestions by Cutri et aL (1985), 
Kinney et aL (1991) and Paltani & Courvoisier (1994), and i) give a quantitative 
estimate of the average behavior of the population as a function of time interval; 
ii) are not biased in favor of highly variable objects, since the objects are extracted 
from the complete flux limited PG quasars. Of course, in this kind of analysis any 
phase information is neglected. 

4. DISCUSSION 

The results of our optical observations can be compared with previous studies (e.g 
Trevese et al. 1994, Hook et al. 1994, Cristiani, Vio & Andreani 1990), which show 
larger variability at all time intervals. In general the variability may depend on the 
intrinsic luminosity, on redshift and wavelength. In particular the dependence on the 
intrinsic luminosity has been recently studied in a large sample by Hook et al. (1994). 
This dependence has also been found by Trevese et al. (1994), although in a smaller 
and intrinsically fainter sample. For this reason, the comparison should be performed 
using objects in a given range of absolute magnitudes. The dependence of variability 
on redshift and frequency is less clear and requires further analysis. In fact, in paper 
I, it has been found a positive correlation of the amplitude of variability with redshift 
and it has been suggested the following explanation in terms of spectral variability. 
A fixed observing frequency corresponds to higher rest-frame frequencies at higher 
redshifts, thus a hardening of the spectrum in the bright phase, i.e. a larger fiux 
variation at higher frequencies, implies a stronger variability at higher redshifts. The 
dependence of variability on redshift has been confirmed by Trevese et al. (1994), 
and is present also when the analysis is restricted to the objects with Mb < —23.5. 
Hook et al. (1994) do not find such correlation in their larger sample. However, 
the fiux limit produces strong correlation between absolute magnitude and redshift. 
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implying that high-redshift objects are intrinsically brighter on average. Thus the 
intrinsic negative correlation of variability with brightness could mask the positive 
correlation with redshift. 

Simultaneous multiwavelength observations of complete QSO samples, with a 
uniform coverage of the luminosity-redshift plane, are required to disentangle the 
various dependencies (see Cristiani et al. 1995). However, also non-simultaneous 
observations of homogeneous samples at different wavelengths can provide new in- 
formation about the behaviour of variability. 

The dependence of the structure function on frequency, which we have found for 
the PG sample (see Fig. 3), allows us to elucidate the dependence of variability on 
frequency and redshift, through a comparison with the optical data from previous 
studies. We have chosen two reference intervals of At, 0.3 it 0.09 yr and 2 it 0.6 yr, 
motivated by the fact that at shorter time intervals the estimate of the variability 
is more affected by the noise, while at longer time intervals high redshift objects 
are progressively lost, due to the finite time baseline. We have then plotted in Fig. 
4 the variability as a function of the average rest frame frequency of each sample. 
The points in Fig. 4 show an increase of variability with the rest frame frequency 
at both time intervals. Other existing variability data are derived from blue band 
samples (Cristiani et al. 1990, Hook et al. 1994, Trevese et al. 1994) which we 
have restricted to 2; < 2.2 and Mb < —23.5 for comparison with our PG data. At 
the relevant redshifts (z) ~ 1.3 — 1.5 the average rest-frame wavelengths lie in the 
interval 1900-1800 A, close to the rest-frame wavelength of the lUF-long camera 
for PG quasars which are on average at lower redshift, (z) ~ 0.3. The average 
variability of the above samples for the two adopted time intervals is shown in Fig. 
4 and appears consistent with the overall trend despite the different distributions of 
absolute magnitudes and redshifts. Finally, we have also included in the analysis the 
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variability data in the [/, Bj and F bands of Cimatti, Zamorani & Marano 1993. 
Due to the hmited time basehne (~ 3 yr in the observer frame) the corresponding 
points are plotted in Fig. 4 only for the shorter rest-frame time interval of 0.3 yr. 
These points appear consistent with the overall trend. In particular the increase 
of variability from the F to the U band is confirmed within the sample, and the U 
point is fairly consistent with the lUF-short data of the PG sample. 

The smaller variability in the red could in part be ascribed to a dilution of the 
quasar light by the non variable emission of the host galaxy. We estimate that this 
effect is negligible in our sample for the following reasons. Indeed Flvis et al. 1994 
evaluated the fractional contribution / = Ig/ hbs of the host galaxy to the observed 
flux. From their Fig. 6 the value of / at log i/^est = 14.26 (A = 1.65/im) appears to 
be ^ 20%, for objects with My < —24 corresponding to our limit Mb < —23.5. We 
extrapolate the correction factor to log Urest = 14.86 using their template spectrum 
for the host galaxy obtaining an average correction / ^ 13%. Since the observed 
flux is lobs = Iq -\- IfJ = Iq -\- fIobs7 then the intrinsic quasar variability 5Iq/Iq = 
(1 — f)~^SIobs/ lobs increases only by a factor 1.15. According to the above template 
spectrum, the correction in the UV is smaller and can be neglected. However, even 
applying the correction only in the red band, the average slope of 5*1 vs log frest 
(Fig. 4) is not appreciably affected. 

The increase of quasar variability with rest-frame frequency implies an increase 
of variability with redshift for samples monitored at a flxed photometric band. 
This effect has been pointed out by Giallongo et al. (Paper I), who invoked it 
to explain the positive correlation of variability with redshift found in a composite 
quasar sample. The slope of a straight line representing data in Fig. 4 is about 
AS'i(At)/A log 1/ ~ 0.28 both for At = 0.3 yr and At = 2 yr. This corresponds 
to the same slope AS'i(At)/A log(l -|- 2;) ~ 0.28 in the variability-redshift relation. 
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assuming that, during the brightness change, the power-law shape of the individual 
QSO spectra log = alog[(l + z)i'obs] + const is preserved. 

In a recent paper, Cristiani et al. (1995) have performed an analysis of the 
variability of a composite sample, consisting of data from Cristiani et al. (1990), 
Hook et al. (1994) and Trevese et al. (1994). The authors confirm the increase of 
variability with absolute magnitude and the positive correlation of variability with 
redshift. They show in their Fig. 5 the structure function in intervals of absolute 
magnitude and redshift. Selecting for comparison the magnitude intervals with 
-26.5 < (Mb) < -25.5, it is possible to estimate A5'i/Alog(l + z) ~ 0.25 - 0.30, 
in good agreement with our estimate from Fig. 4. Thus, as suggested in Paper I, 
the entire redshift dependence of variability can be explained in terms of spectral 
variability. 

Recently, gravitational microlensing from a population of intergalactic objects 
has been studied as a possible cause of quasar variability (Hawkins 1993, 1995, 
Schneider 1993, Alexander 1995). Achromaticity is a characteristic signature of 
microlensing of stellar images in the case of "local" galactic lenses. However, in the 
case of intergalactic microlensing the Finstein ring of the lens can be smaller than the 
source size. In the standard accretion disk model the high and low energy photons 
are emitted from the inner and outer parts of the disk, respectively, and this results 
in a color change during the amplification. While it is difficult to distiguish between 
this phenomenon and the intrinsic color variations (Baganoff & Malkan 1995), a 
statistical analysis of the variability as a function of wavelength can constrain models 
of microlensing-induced variability. 

Concerning the intrinsic dependence of variability on frequency, Paltani & Cour- 
voisier (1994) give an explanation in the framework of thermal emission models. 
Assuming that the variability is due to temperature changes, a hardening of the 
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spectrum in the brighter and hotter phase is expected, provided that the turnover 
in the spectral energy distribution is not at much higher frequency than the observ- 
ing band. If the temperature is too high, the slope of the spectrum in the observed 
optical/UV band becomes less dependent on temperature thus less variable. The 
increase of variability with frequency which we have shown to be present from the 
red to the UV can constrain the maximum temperature of any thermal emission 
model. 

We remark that if brighter objects are hotter on average and the variability is due 
to temperature changes they should be less variable. Indeed, assuming for example a 
spectrum of black body or thermal bremsstrahlung, the spectral turnover of brigther 
objects is progressively shifted at higher frequencies producing progressively smaller 
flux changes from 5I/I oc 5T /T[hi'/kT) for kT < hu, to 81 j I oc 5T/T for kT > 
hu. This could explain the observed negative correlation of the variability with 
luminosity. 

5. SUMMARY 
We can summarize our results as follows: 

1. The variability of 30 PG quasars has been observed in the red band during 
three years, with typical accuracy of 0.03 mag (Fig. 1), using a Schmidt with 
a CCD detector. 

2. The analysis of the structure function in the rest frame shows that the vari- 
ability increases from the noise level at a time interval At ~ 1 month to ~ 0.2 
mag for At ~ 2 yr. 

3. Using public lUE data of a sample of 21 PG quasars, with 12 objects in common 
with our sample, we have compared the variability in the red and in two UV 
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bands as a function of time interval. The variability appears to increase with 
the rest-frame frequency at each time interval. 

4. The average variability of other published quasar samples has also been com- 
puted as a function of their average rest-frame frequency. The results support 
the evidence for an increase of variability with frequency. 

5. The positive correlation of variability with redshift found by Giallongo et al. 
(Paper I) and recently confirmed by Cristiani et al. 1995 can be entirely ac- 
counted for by frequency dependence as already suggested in Paper I. 

6. Both the correlations of variability with frequency (and hence redshift) and 
intrinsic luminosity are qualitatively consistent with thermal spectra becoming 
hotter, and locally harder, in the bright phase. 
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CNR and MURST. 
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FIGURE LEGENDS 

1 Distribution of photometric errors. 

2 Individual magnitude variations \5m\ as a function of the rest frame time in- 
ter vaL 

3 Structure functions of the PG objects studied in this work. GIMP: Campo 
Imperatore R band; lUE-longA: 2500-3000 A; lUE-shortA: 1300-1900 A. The 
abscissa of each point is the average time interval of the data in the relevant 
bin. 

4 Variability as a function of the rest-frame frequency, for two different time 
intervals. White symbols: 0.3 yr; filled symbols: 2 yr. 
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